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Monte-Carlo sampling of self-energy matrices within sigma-models 
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Abstract 

The 'Neumann-Ulam' Monte-Carlo sampling is described for the calculation of a matrix inversion or a Green 
function in case of Hubbard-Stratonovich (HS-)transformed coherent state path integrals. We illustrate how 
to circumvent direct numerical inversion of a matrix to its Green function by taking random walks of suitably 
chosen matrices within a path integral of even- and complex-valued self-energy matrices. The application 
of a random walk sampling is given by the possible separation of the total matrix, e.g. that matrix which 
determines the Green function from its inversion, into a part of unity minus (or plus) a matrix which only 
contains eigenvalues with absolute value smaller than one. This allows to expand the prevailing Green 
function around the unit matrix in a Taylor expansion with a separated, special matrix of sufficiently small 
eigenvalues. The presented sampling method is particularly appropriate around the saddle point solution of 
the self-energy in a sigma model by using random number generators. It is also capable for random sampling of 
HS-transformed path integrals from fermionic fields which interact through gauge invariant bosons according 
to Yang-Mills theories. 
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1 INTRODUCTION 



1 Introduction 

1.1 Coherent state path integrals and the HS-transformation to self-energies 

Coherent state path integrals with a quartic interaction of fields, especially those with a fermionic character, 
allow for the transformation to even- and complex- valued, anomalous-doubled self-energy matrices PQ-[I]. After 
such a Hubbard-Stratonovich (HS-)transformation [6], one usually performs a coset decomposition of the 
total self-energy in combination with a saddle point approximation [7]- [9]. In most cases the self-energy matrix is 
restricted to the elements of the orthogonal so(9T, 91) or symplectic sp(9t, 9t) Lie algebras for fermionic or bosonic 
degrees of freedom in the initial coherent state path integral, respectively. Therefore, the transformed, total self- 
energy matrix can be factorized into density-related and anomalous-doubled parts with a coset decomposition 
80(01,01) /U(9t) ® U(9t) for fermions or Sp(9t,9t) / U(9t) <g) U(9T) for bosons. This results into a gradient 
expansion of the determinant with the coset matrices which precisely contain the 'Nambu'-doubled, coset 
matrix algebra field degrees of freedom for a spontaneous symmetry breaking |10[ [TT] . 

In this article we point to the remarkable property of coherent state path integrals concerning the 'Neumann- 
Ulam' Monte-Carlo sampling of Green functions |12j-|14j. The presented case avoids direct numerical compu- 
tation of a matrix to its Green function, as e.g. by recursive techniques for approximate inversion of partitioned 
sub-matrices (cf. Refs.[T5]-|18]). We describe the applicability of random walks within suitably chosen matri- 
ces occurring in coherent state path integrals so that matrix inversion is obtained from random sampling of 
transition probabilities of matrix elements. This is combined with properly defined scores of residual weights. 
The chosen system of Refs.pjJJ [20] (section 3.3)H is generic for other path integrals with self-energies as that 
of QCD so that lattice QCD sampling with determinants can be replaced by the presented 'Neumann-Ulam' 
Monte-Carlo method which is also applicable to the HS-transformed coherent state path integral of BCS-QCD 
in Ref. |21l 122^1 . By way of example we consider the model of Refs. [T9l [20] (section 3.3) with a two-band, 
semiconductor-related solid in order to hint to the possible 'Neumann-Ulam' sampling. The corresponding 
HS-transformation of the original fermionic coherent state path integral is achieved from dyadic products of 
anomalous-doubled, anti-commuting fields which results into an even- and complex-valued self-energy matrix 
as an element of the so (91, 91) Lie algebra. The presented sampling method is particularly applicable around the 
saddle point solutions of the anomalous-doubled self-energy so that the fluctuations can be determined around 
the saddle point approximation with random number generators. The implementation of a random sampling 
for HS-transformed coherent state path integrals follows from the possible splitting of the total matrix, e.g. 
that matrix which comprises the Green functions from its inversion, into a part of unity minus(or plus) a matrix 
which only consists of eigenvalues with absolute value smaller than one. In consequence, the prevailing Green 
function can be expanded around unity in a Taylor expansion with the separated matrix whose eigenvalues are 
sufficiently confined within the open interval between minus and plus one. 

1.2 'Neumann-Ulam' matrix inversion 

In the following we briefly review the 'Neumann-Ulam' matrix inversion in order to extend this sampling method 
to HS-transformed coherent state path integrals, as e.g. given in Refs.[l9l 120] and |21[ 122]. The 'Neumann- 
Ulam' Monte-Carlo sampling takes random walks of a matrix Mij in order to achieve a matrix inversion M^ 1 , 

as for a Green function G%j = Mrr . However, this method can only be applied, provided that the matrix 
can be separated into a unit matrix ljj minus a matrix rhij whose eigenvalues are completely restricted to the 

2 Title: 'Fermionic coherent state path integral for ultra-short laser pulses and transformation to a field theory of coset matrices' 
3 Title: 'BCS-like action and Lagrangian from the gradient expansion of the determinant of Fermi fields in QCD-type, non-Abelian 
gauge theories with chiral anomalies' 
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open interval ( — 1, +1) so that following expansion holds 

Mij = lij — rhij ; (eigenvalues of rhij E (—1, +1) ) ; (1.1a) 
Gij = M^ 1 = (i - m)". 1 = lij + rhij + rhii 2 rh i2j + rhu 2 rh i2 i 3 rh i3 j + . . . . (1.1b) 



The 'Neumann-Ulam' method uses a decomposition of the matrix rhij = Pij Wij into transition probabilities 
Pij > and a fixed, residual weight Wjj G C so that random walks within the matrix rhij = Pij of various 
order (rh k )ij result into the matrix inversion of M^ 1 

Gij = M^ 1 = (i - mjr. 1 = lij + (pij wy) + (p ii2 <N ii2 p i2j <N i2j ) + (p ii2 <N ii2 p i2h <N i2is p hj w i3 j) + ... ; (1.2a) 
rhij = Pij . (1.2b) 

In order to obtain a valid score, one has also to introduce a stopping probability Pj = 1 — Y2j Pij so that 
a particular random walk does not continue 'ad infinitum' within above series (jl.lbp and (jl.2ap ! A valid 
random walk {(rh k )ij) of order k within above summation of Gij = has therefore to include the stopping 
probability Pj into the 'total probability density J" of transition probabilities and an additional factor 1/Pj into 
the 'score S' of the residual weights so that the decomposition of rhij = Pij still holds for a valid expansion 

"i^/,, 1 (i '"),/ 

((™%-) = ( (p«2 Pi 2 * 3 ■ ■ ■ Pui fc+ i p i fe +i) X (*« 2 w; 2 ; 3 ■ ■ ■ < N i k i k+1 b— K+id) ) ; ( L3a ) 

Xv * ' Hfc+l ' 

'probability density' IP v v ^ 

'score' §| (™«ini]r) 

P i = l-^ Pii . (1.3b) 

As we compare the well-known sampling of an integral / = dx /(x) with a normalized probability density 
p(x) > and residual weight function w(x) = f(x)/p(x) to the presented case of random walks with a matrix 

rhij 

+oo r+oo 

dx f(x) = j dx p(x) w(x) ; (1.4a) 

-oo J — oo 

+oo 

dx p(x) ; p(x) > ; /(x) = p{x) w(x) , (l-4b) 

-OO 

one is guided to the proper definition of a probability measure in analogy to 'dx p(x)' and to the definition 
of the corresponding score for every random walk of order k. Note that the overall, total probability measure 
CP{p, P} of random walks of all lengths is given by the sum of its sub-orders of definite length m. Every sampled 



sub-order J , i™ ) im+1 {pr S , P m+ i} oc p hi2 p i2 i 3 ■ ■ ■ p imim+1 ('m') selects a particular sequence vj ili2 \% 



«2*3 ' ' ' W im'ra+1 



which is abbreviated by the specific score (w rs ; ^-) with the same length J k = m' and the same identical 
path of residual weights in between the start- and end-points i\ = i, i m +± = j. Since a random walk cannot 
continue for infinity, one has to insert the corrections of the stopping probability P« m+1 and 1/Pj into probability 
measure and score, respectively 



7{p, p} = £ y£l im+1 {Prs, P, m+ J ; ?Sr 0) {p-, p ul = i ; (i.5a) 

m=0 
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^ > il...i m + l{P»".S! Rlm+1 } — Plll2 P«2«3 ''■ Pimim + l Pim + 1 > (1.5b) 

assigned score : => S^(w rs ; 4-) = 6i lt i <Ni x i 2 Wj 2 j 3 ■ ■ ■ Wj fe j fe+1 h k+1 ,j &ij °\^rs', p") = Uj ■ (1.5c) 

The 'Neumann-Ulam' averaging for the Green function Gij, abbreviated by following brackets {...) N .u. , is 
therefore properly determined by using the probability measure ifc+i {p rs , Pj fe+1 } and its accompanying 
score sj^(w rs ; p-) for every order T of the expansion (1 — f^)^- 1 

oo 

Gij = Mi/ = (i - m)^ = ((i - w)" 1 ) = l {j + £>{p, P} §J'(w rs ; = (1.6) 

' k=l 

oo OO ^ 

= + 2 ^Swi^"' p ™+ii §!j fc) (Wrs; pr) = ly + J] (P*a Pfcis ' ' ' Ptmj P j) («% Wtafa ' ' ' ^mj p 

k,m=l m=l J 

Hence, the variance straightforwardly follows by application of the probability measure fll.5all.5bh and the score 
(|1.5c[) and by introduction of an additional matrix K{j = rhij Wij which determines the deviation from the mean 

value combined with the inverse of the stopping probability 1/Pj. We emphasize that each score S^(w rs ; 4-), 

S^(w rs ; pr) of order 'fc' and T is specified by a sub-probability measure 7^ ik+i {p rs , Pj fc+1 }, 7^ i[+i {p rs , P« i+1 }, 
respectively. As one takes a random walk of length 'm' according to the transition probabilities and stopping 
probability with the sub-probability measure im+ii-Prs, Pi m+ i}) one is a l so naturally confined to products 
of scores S,-^(w rs ; 4-), S,-?(w rs ;4-) which have the same length W and the same identical path of transi- 

"J y 3 J 3 

tions of the actual, sampled random walk with 3\-™\ m+1 {p r s, Pj m+ i}- in consequence the sub-probability mea- 
sure J > ,-™\ m+1 {prs, Pj m+ i} always specifies the actual random walk and restricts the allowed, possible scores 
S^(w rs ; pr), §>fhw rs ; 4-) to the same identical steps and the length of the prevailing sub-probability measure. 

" 3 •* 3 

We have especially focused to this point because this leads to the important, two Kronecker-deltas 5 m ^ and 5 m i 

in following expansion for the variance with the sampled path ^™.i m+1 {Prs, P m+ i} = Phi 2 Pi 2 i 3 ' ' ' PimWi P Wi 
including the correction of the stopping probability 

(AG,) 2 = (((i- w) 5 ' - «i - *tfW)'L. = (<i - *tf (i - - (((i - *tf>J' - 



£ ((«*)(,■ (*■)«) -(G„) 2 = Sjj'l*"^' S«(w„; ^) - (G„) 2 = (1.7a) 

fe,«=0 ' ' fc,Z=0 

OO OO 

' P iI n) .im+i^P rs ' P Wli ^m,fe S!f(Wr S ; 4") ^m,i (w rs ; p-) - (G^) 2 = 



fc,i=l m=l 

oo \ 2 

+ ^2 3> ff< m+ l'fP r< " Pim + ^ ^ ^ 2 ' ( Wi2 * 3 ' ' ' ^Wl) ' ( p ) ^m + l,i - (Gij) 2 

m=l 

lij + ^2 ( Pli2 P* 2 *3 ' ' ' P*^' P J') (™«2) 2 (Wj 2 j 3 ) 2 • • • (Wi m j) 2 ( p") - (Gij) 

m=l ^ 

°° 1 1 

Uj + ( £ k ^ k ^ ■ ■ ■ Kimj) p. ~ (G*jf = (1 - pT - (^) 2 

m=l ^ ^ 
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Kij = pij (w/jj) 2 = rhij <Nij . (1.7b) 

A variant of this sampling method is given by a partial removal of the stopping probability Pj. If the random 
walk does not stop after 'k' steps according to the conjugate probability Pj = 1 — Pj = ^ • for continuation, 
one simply adds the actual, obtained score w/jj 2 Wj 2 j 3 • • • Wj fe ,- without the inverted stopping probability to the 
already existing, accumulated value of order 'A;' and generates another transition probability Pi k+1 i k+2 for a 
further step 'k + 1'. One can thus proceed to higher order scores '(w')', (/ > /c) with the actual random walk 
which only ceases as soon as a stopping probability P/ is sampled in a higher order 7 > k 1 of the random walk. 
The corresponding probability measure and score has to be modified in order to take into account the partial 
removal of the stopping probability Pj with the conjugate continuation probability Pj. We therefore distinguish 
between the two cases Qj = Pj and Qj = Pj for stopping and continuing of the random walk in following modified 
definition of probability measure and score of a random walk 

Pi = i - fa ; Pi = i-Pi = J2 fa ; ( L8a ) 

3 3 

oo 

?{p,p,p} = i(np,P} + np,P}) = \ J2JT 9^ dm+i {p rs ,Q im+1 } ■ y^itriM = i ; 

q =p pm=0 

(1.8b) 

{P rs ' Q*m + 1 } = Pn«2 P«2«3 ''' PWm+1 Qim + 1 ! (Q«m+1 = Pim + l OT Q«m + 1 = Pim + l ) > (1.8c) 

§i?(wr S ; qt) = Wi!i 2 Wj 2 , 3 • • • w, fcifc+1 -— S ih+ld ; S^ =0) (w rs ; ±) = Uj . (1.8d) 

According to relations (jl.5alll.6"j) . one can define the modified 'Neumann-Ulam' sampling by the brackets 
(. . .)^ with the probability measure ?{p, P, P} = ±(3>{p, P} + ?{p, P}) ()1.8aJ1.8cj> and score §if{w rs ; ^-) 

oo 

G l3 = Mr 1 = (i - m)r; = ((1 - w)- 1 ) _ = + \ *{P, Ql s S?(wr S ; = (1.9) 

Q = ppfc=l 

OO 

Q=P,P fc,m=l 

oo ^ 

= + 5 E E (p" 2 pi2i 3 ' ' ' PW Qj) (^«a ^2« 3 • • • Wj m j q- 

Q = ppm=l ■? 



lij + 5 E [(P«2 P ^3 ' ' ' P*mj P j) (w« 2 Wj 2 j 3 • ' • Wj m j — ) + 
m=l ^ 

+ \Pii2 P«2«3 ' ' ' Pimj (1 _ Pj 



^ (1 - Pj). 

In a similar manner we can conclude for the corresponding variance of the modified random walk (. . .)fpu with 
a partial removal of the stopping probability. As one applies the proper definitions (ll.8alll.9j) for the sampling, 
one achieves a variance with an additional term which is proportional to the inverted continuation probability 
1/Pj = 1/(1 - Pj) (cf. (jl.7all.7bj> ) 

(AG,/ = (((i -w) 5 ' - (<!-*#) ) 2 ) = ((i-w) 5 ' (t-*)-'\ - «(i -*)«') J' = 
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= I 1 - \Qr + (Tr^y ) - ^') 2 ; k * = p*j 0%) 2 = ^ W y • ( L1 °) 

One might object that the possible split Mij = lij — rhij already simplifies the calculation of the corresponding 
Green function Gij = M^ 1 = (1 — rh)^ 1 = lij + fha 2 rhi 2 j + ... in such a manner that it is sufficient to compute 
matrix multiplications (fh k )ij at varios powers 'k' without any Monte-Carlo sampling or random walks of 
transition probabilities. This may be advantageous for strictly banded, 1-dim. cases; however, as more and 
more matrix multiplications are performed to higher orders (fh k )ij, the spread of an initially sparse or banded 
matrix rhij immediately leads to rather occupied, full matrices which tend to the full amount of multiplications 
iVj, (rhij]i,j = 1, . . . ,Nrh) at every order 'k' in the Taylor expansion. Therefore, the matrix multiplications of 
the Taylor expansion up to maximal power 'A; ma:r ' without random sampling will finally approach the number 
'~ k max iV? ' of floating point operations even for an initially sparse matrix. It may then even be more adequate 
to determine the full matrix inversion Gij = (1 — rh)^ 1 without any Taylor expansion by 'direct' numerical 
inversion which also is of the order of '~ TV? ' floating point operations. Hence, we can point out the advantage 
of the presented Monte-Carlo sampling for 'Neumann-Ulam' matrix inversion of Green functions; the advantage 
of a vanishing spread of the initial (sparsely distributed) matrix entries of rhij is kept throughout all random 
walks as higher orders of the Taylor expansion are considered up to the end of the individual random walk with 
stopping probability Pj. The reliability and effectiveness of the described Monte-Carlo sampling for a matrix 
inversion finally depends on the statistical error which is determined by the matrix Kij = rhij and the 
stopping and continuation probabilities Pj, P, = 1 — Pj in the variances (jl.7a|1.7b|) and (ll.lOp . respectively. 

2 Random walks applied to HS-transformed coherent state path integrals 

2.1 Self-energy S^ S1 (T^ 11 ,^;^ 32 ,^) of densities 'a = V and of 'Nambu'-parts 'a ^ V 

In Ref . [19] we have outlined the derivation of the path integral in terms of the anomalous-doubled self-energy 
which one attains after a HS-transformation of the dyadic product of 'Nambu '-doubled, anti-commuting fields. 
In this section we start out from the final result of section 3.3 in Ref.|19j and introduce additional abbreviations 
for the contour time steps T, 3 , the spatial coordinate vector x, the electron-, hole- (/i = e, h) and spin-labels 
(s =t, i) of the semiconductor-related solid. Furthermore, we have to include 'Nambu'-indices 'a, b, c, . . . = 1, 2' 
for the anomalous-doubling of the fermionic fields where the block-diagonal parts '11', '22' and off-diagonal parts 
'12', '21' of the self-energy matrix E"^ Sl saOT,/ 1 ,xi]7 j2 is ,£2) specify the density-related and anomalous- 
doubled parts, respectively. We extremely simplify notations in order to display the structure of the path 
integral with the self-energy matrix as an element of the so(9T, 91) Lie algebra. The overall dimension '9T 
of latter Lie algebra so(91, 91) is determined by the product of the respective two degrees of the band-label 
({j, = e,h), the spin-label (s =t>4-) as well as the time contour index (rjj = ±) multiplied by the number of 
chosen discrete time steps (N = t^/At) and the number (Af x = (L/Ax) d ) of spatial points on an underlying, 
'd'-dimensional grid of system length 'L' with space intervals Ax 

dimension 71 of so(%Vl) = (m = 'e\'h') x (s =t,|) x {m = ±) x (N = t N /At) x (N x = (L/Ax) d ) . (2.1) 

V v ' V v ' v ' V v ' V v ' 

2 2 2 N A4 

We abbreviate the set of variables 'contour time 7 3 ■ 1 ', 'spatial coordinate vector x\ 'band-label (fi = e,h)' and 
spin index '(s =1\4)' by the uppercase letter '3' and the restriction to 'contour time 7 j 3 ' and 'spatial vector 
X* by the lowercase letter i 

(7^]x^,s)=3; (T^|x) = i. (2.2) 
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The notation of the self-energy matrix S1 ^(T^ 31 ,xi;T^ 12 , X2) is thus transformed to E^ 6 



J 2 



^A<i.aiiA»2,*a(*^x J1 >^l! ^2 IS '^2) -> S ai3 2 ' ( 2 - 3 ) 

where additionally attached subscripts '1', '2' of '3i', 'J2' (or respectively the reduced case Hi', H2' ) distinguish 
between the two sets of contour-time-coordinate labels including band- and spin indices for the uppercase 
letters (|2.4a|) or respectively lowercase letters (12.4bp 

3i = ; ^2 = (T^ 2 !x 2 ;/i 2 ,S2) ; ■■■ ; (2.4a) 

i 1 = (7^]x 1 );i 2 = (7^]x 2 ); .... (2.4b) 

According to the dyadic product of anomalous-doubled, fermionic fields, the orthogonal symmetry 80(01, 9T) 
causes 'Nambu'-parts (b ^ a) |23[ [24"] of an overall, hermitian self-energy matrix E^ whose diagonal block 
parts are related by opposite sign and transposition and whose two, anti-symmetric 'Nambu'-parts are related 
by hermitian conjugation 

/ yll yl2 \ab 
yab _ ^Cfida ^lSa \ ■ yaa _(yaa NT. y22 __(yll \ T . 

^3132" I £21 y22 J ' ^OWa - V^iO^ ' ^3i?2~ K^xJn) ' 

^3 2 =(^a 2 ) 2t ; ^=-(sg 2 b) ) T ; tf i3 =0 (for( J1 wj 2 = 0)^(^^2 = 2^ + 1)). 

(2.5) 



A term as 



V Tr 

^ a,b 



(E^/AQ ■ (S^/AQ , (2.6) 



with the distinguishing subscripts '1', '2' of '3i', 'J2' (|2.4a[) therefore implies a summation of the contour 

(n ■ ) (77 ■ ) 

times 'T,- Jl ', "J" J2 J2 ', the spatial coordinate vectors 'x\\ 'X2', the band- and spin-labels , Hi, 112', ' 81,82' apart 
from the supplementary trace over 'Nambu'-space (a, b = 1,2) with the anomalous-doubled components of the 
self-energy. Note, that we use normalized, spatial summations with factor 1/J\f x which we can absorb by a 
rescaling of the self-energy within our derived path integral of section 3.3 in Ref.|19|. (The parameter M x of the 
total number of spatial points is similar to the parameter of the matrix dimensions in random matrix theories 

(SUES].) 

tf l3a /Af x -+ ±f l32 . (2.7) 

Aside from the sets 'J2', Hi', H2', (|2.4al2.4bp of collected variables, we abbreviate the charge 

index 'q^ = ±1', 'q^ 2 = ±1' of electrons (/ii,^2 = e) and holes {^1,^2 = h) and the contour time index 
'rjj 1 = ±', 'rjj 2 = ±' in a similar manner to relations (I2.2l2.3p by subscripts '1', '2' 

gi = <?m ; 92 = ; • • • ; ( 2 -8a) 
m = Vh ; V2 = vh ; ■■■ , (2.8b) 

in the overall summations of the contour-time and space coordinates with the band- and spin-degrees of freedom. 
Moreover, one has to regard the 'Nambu'-metric tensors I afe , S ab which have to be inserted for the proper dyadic 
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product of fermionic fields as the HS-transformation is performed for the remaining self-energy of the 

so(9t, 91) Lie algebra 



I ab = 5 ab diag{ 1 , i } ; S ab = 5 ab diag{ 1 , -1 } = l ac ■ l cb ; (a, b, c = 1, 2) . 



(2.9) 



0=1 



a=2 



a=l 



a=2 



After a rescaling of the one-particle operator H§j an d the source term 3^j 2 as for the self-energy EJJ* 



a6 



^06 



J 2 



ab 



Tab 



nab 



(2.10) 



the ensemble-averaged generating Z[3] of section 3.3 in Ref.flU] can be specified in terms of a notation which 
signifies a similarity to random matrix theories 



Z[3] 



x DETi 



d\£ ba 



a' ,'2' 



x Tr 

a,b 



(2.11) 



At 



ba 



+ 7 



Sections 2.1 to 3.3 in Ref.|19] consist of the derivation with the HS-transformation from a path integral of 
anti-commuting fields for a semiconductor-related solid coupled to an ultra-short laser pulse with an inter-band 
dipole moment. The latter, driving light-semiconductor interaction causes off-diagonal '12', '21' 'Nambu'-blocks 
in the one-particle operator Hj apart from the usual density-related kinetic energy terms with further 

possible one-particle potentials along the block-diagonal parts '11', '22'. The Coulomb interaction V(xi,X2) 
and the correlation function f(x\, ii; £2, £2) of disorder and noise are combined into the generalized interaction 



potential V i2il =V(x 2 ,T^ 2) ;x 1 ,T;; ,j1 



1 =1 

2 h 



ul /(f 2 ,'J 72 ;fi,'J il ); (ji,j 2 



1, 



,2iV) 



(2.12) 



which enters by the quotient (V^) 1 for the quartic interaction of Fermi fields into the Gaussian factor with 
the two self-energies fij^j , ^j^ 2 • The latter matrices are accompanied by the charge- 'gi, 52' an d contour-time 
'771, 772' indices including the diagonal 'Nambu'-metric tensors S aa , S bb . As we apply the simplified notations 
(I2.21l2.8b"|) and the rescaling of the self-energy matrix Sj j (12. 7p . of the the one-particle operator Hj"^ an d 
of the source term 3y 2 3i lOj) to the derived path integral (|2.1ip from section 3.3 of Ref.|19j , one obtains 

following expression Z[3] f)2. 13j) of the generating function whose 'exponential trace- log' relation replaces the 
determinant in fl2. 1 If) . The pre-factor \ is caused by the square-root operation of the determinant in (12. lip 
with the additional sum over the discrete contour time steps beginning at j = and ending at j = 2N + 1. 
We have separated the total number of spatial points M x as a parameter for the saddle point approximation 
by using a normalized spatial sum over the discrete grid points 



Z[3}= / d[£ 



%r x ] exp 



a, 6 - 



(2.13) 



x cxp 



E- 



2JV+1 

iAA, V V Tr In 



a,6 



bb Ljfca 



I H 



:, aJa I aa + %a 32 3^i + S fe6 92% 



At 



x sum over the discrete spatial grid points 
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The familiar kind of the generating function Z[3] (|2. 1112.1312, 14a|) for a sigma model gives rise to a saddle point 
equation (|2.14b|) with solution si, ^ a from a first order variation aS^ 3i ()2. 14c|) with respect to the original 
self-energy Sj" 3 



Z[3] = I d[X b £ 3l ] exp 
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(2.14a) 

(2.14b) 
(2.14c) 



We disregard the detailed model and theory of Refs.[19j-[22j with the semiconductor-related solid coupled to 
ultra-short laser- light pulses and perform an expansion of the generic sigma- model Z[3] (|2. 1112.1312. 14ap around 
its saddle point solution Sj°^ a of the self-energy by using the corresponding Green function (|2.14b|2,14cp . 

As a result we can sample random numbers for AE^j around the mean values S^^ a for fluctuation properties. 
The kind of saddle point approximation and generating function is generic for other kinds of coherent state 
path integrals as lattice QCD. After the expansion around the saddle point solution Sj°j , we finally acquire 

the generating function Z[3] (|2,15all2.15dp in a manner which, according to the description of the 'Neumann- 



Ulam' Monte-Carlo sampling in section [l~2l contains the crucial decomposition Gij = = (1 — into 
unity minus (or plus) a part 'm oc (\At/h) AE^j ' (|2.15dp whose eigenvalues tend to zero for sufficiently 
small time steps 'oc (\At/h)\ The presence of the saddle point solution in 'rhj^ (l2.15dp . which multiplies 
the fluctuation term 'oc (\At/H) aEj ^ ' in the determinants, does not alter this scaling 'oc (i At/fi)' to smaller 

values because the saddle point solution S^^ a as a mean value approaches a fixed, constant value for physically 
smaller sizes of 'oc (\At/H)\ This property can be strictly proven in an analytic manner by investigation of 

the saddle point equation (|2.14bp so that the terms as 'rhj^ = Q^J^ $ bb 92 V2 ('if) A ^'3' C J3 ^ 3 Q2.15dp 
also advance towards zero values for smaller time steps 'oc (\At/h)\ In the following tabulated relations 
(|2,15aM2.15dp . we point out the similarity to the described 'Neumann-Ulam' sampling in section 11.21 with 
the important split M = 1 — rh (| 1 . 1 atil . 2b|) by under-bracing the determinants DET, their corresponding 
'exponential-trace-log'-relations' DW£(matrix) = exp{Tr ln(matrix)} 7 and the inversion of 'M' in terms of 

>rh jlj3 = §53f S 6b 92 m (i^)AS^ 33 S cc 93 r/ 3 ' (|2~T5djl . Note, that the linear term of 'm^,' (p~~5dl) has to be 
cancelled in the expansion of the logarithm, due to the saddle point approximation (j2.14ati2.14cp which sets 

linear terms of aS^ in the exponents of the derived generating functions Z[3] (I2.15all2~.15bp to zero. This 
latter peculiarity is emphasized by a tilde above the determinants and the logarithms within the 'exponential- 
trace-log'-relations' DET (matrix) = exp{Tr \n(rnatrix)y (cf. 12. 15c)) 



Generic case around mean value j" with sampling of aSj" 3i 
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lll.lall.2bl l: [1- 



linear 'AS^jj '-term has to be cancelled in the exponent of 'DET (matrix) = exp{ Tr \n(matrix)Y I j 



(2.15c) 
(2.15d) 



The condition on m^- (jl.laHl.2b"]) in Gij = M^ 1 = (1 — rh)^, and the restriction of the eigenvalues to be within 
the open interval (— 1,+1), is particularly applicable for the HS-transformed, coherent state path integral 

Z[3] (j2.15aM2.15dp because the analogous self-energy matrix or term in place of m^- (j2,15dj) turns out to be 
proportional to the time step \At/H times the self-energy matrix ASfe^ = AS^ S2;/iijSi (t] 2 J2 | x\) itself 

with incorporation of the Green function factor §3^3 ()2.14bj) . The condition on the eigenvalues of rh cx 
(\At/H) AS is thus always achievable for sufficiently small, time-step intervals At which has anyhow to be the 
condition for proper, analytical and numerical considerations. As one is inclined towards smaller time steps with 
the analogous self-energy terms rh oc (i At/H) AS for a fixed time to, t\,... , t/v, (At = tjy/N) of observation, one 
attains smaller eigenvalues of rh oc (i At/H) AS and a better convergence of the expansion M^ 1 = (i — m)" 1 , but 

has to carry out the random sampling for higher dimensional matrices aSJ"^ = Al& S271 S1 (T^ 2 ] X2; T^ 1 ] x\). 
However, larger time steps 'At' reduce the effort of the sampling of the random, 'fluctuation' self-energy for a 
fixed observation time, but may corrupt the condition on the eigenvalues of rh oc (\At/H) AS to be within a 
sufficiently small interval of (— 1, +1). Too large values of 'At' can even result in a total failure of the general 
time development of the considered model. Therefore, a crucial choice in between the two described extremes of 
rh oc (\At/H) AS may lead to the most adequate sampling around the saddle point for a practicable numerical 
investigation. The condition on a proper value of 'At' has to be combined with the expression of the variance 
(ll.7all.10j) which includes the matrix Ki 



rrii 



Wij (jl.7bll.lTJl) . Larger transition probabilities Pij imply smaller 
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residual weights and hence a smaller value of (1 — K)^ , but increase the inverted stopping probability 
1/Pj, resulting to an overall larger variance. Smaller transition probabilities pij generally cause larger residual 
weights Wij and therefore a larger value of (1 — K)^ , but decrease the inverted stopping probability 1/Pj. In 
consequence a good result of the 'Neumann-Ulam' sampling for (|2.15ati2.15cp requires a subtle compromise in 
the choice rhij = Pij Wjj of the transition probabilities and of the scores with the residual weights in order to 
achieve a small variance (ll.7ap . Similar considerations hold for the modified 'Neumann-Ulam' sampling with a 
partial removal of the stopping probability (jl.lOp . 



3 Conclusions 



3.1 Choice of probability densities and approximation schemes of the random sampling 



An extension of the 'Neumann-Ulam' matrix inversion can be proposed for rather general functions [/(M)]j 



[/(I 



rh)]ij of a matrix Mij 



h 



m, 



, provided that the eigenvalues of rhij are sufficiently placed near zero 



so that a Taylor expansion of [f(l — m)] converges with respect to the matrix rhij. Neglecting the first line of 
Z[3] (I2.15ati2.15bp with the exponential of the actions from the saddle point solution T>~°\ ab , we can separate 



^ba 



the remaining integrand in the path integrations over AS^ 3i (j3.lati3.lbp into the probability density $H[aE^ 2Ji j 
(|3,lcp . accompanied by a phase factor ^[aE^jJ (|3.1dp . times the Green function &j 3 j 1 (|3.1ej) which can be 
computed according to the 'Neumann-Ulam' matrix inversion because of the given split into (1 — m)^ 1 with 

™hrs = ^y^2 ^ ^ 2 Of) A ^3233 ^ 3 ^ 3 Q2.15d)3.1f|) . In consequence one has to take random numbers 
corresponding to the distribution (ifAS^J D^AS^aJ which, itself, can also be tremendously simplified and 

approximated according to the split '1 — m' in the determinant DET( [1 — m] ) + 5 beginning with a Gaussian 
factor exp{— ^Tr[ (m) 2 ] (cf. (|2.15c|3.Tg| ))! It is possible to shift remaining higher order terms of traces to the 
Green function (S^^ (j3. le[) which yields a modified observable. Then the choice of transition probabilities and 
residual weights has to be determined so that one can sample random walks for the inverted matrix or Green 
function &"^ 3l following from the described 'Neumann-Ulam' matrix inversion method of section 11.21 
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3 CONCLUSIONS 



=m[l- 9$? S bb a 2 m (if) a£^ 3 S- , 3 % ] ^ sg>f * = * [l - rh) §gf * ; (3.1e) 
Sua = SSf ^ Q2 m (if) Alfe, S cc q 3 m ; (3.1f) 
DET( [i - m] )+* = exp { ^Tr ln[i - fh] } = exp {- i Qtt[ (m) 2 ] + ^Tr[ (m) 3 ] + ^Tr[ (m) 4 ] + •••)} (3.1g) 

A possible choice of transition probabilities and residual weights Wjj for a sparse matrix = p^j Wj» ()3. If [) 
should result from the factorization 



Pij = s J\rhij\; \wij\ = J\rhij\/s ; rhij/\ihij\ = Wii/|wy|; my, e C , (3.2) 

with an appropriate scale factor sq so that a direct numerical inversion can be circumvented. A good choice of 
stopping and continuation probabilities Pj , Pj = 1 — Pj should be attained for equal values of | from the relation 
Pj + Yli Pij = Pj + s o Yli vPhTI = 1 with rhij (|3.1f j) being proportional to the time interval At and to the 
saddle point fluctuation aE^ . A final adaptation with the overall scale factor sq determines the transition 
probabilities and residual weights. A further choice of transition probabilities can include the physical intuition 
where higher precedence is given to those states "3x\ '^2' of the transition probabilities 'fb^' an d residual 
weights 'w3 1 3 2 ' which yield a higher occupation probability of particles in the physically important parts of the 
spacetime coordinates. 

Nevertheless, there is a further problem of the random sampling according to the probability density 
rf[AE§*jJ £H[a£^J (IXTc]) which, however, occurs in any Monte-Carlo method with 'many' independent ran- 
dom variables. If one reduces the probability density ^[aE^J ^IA^mJ P-1 C P to the analogous problem of 
'N z ' independent, random Gaussian factors 

N z 

d[Aflg ai ] ^AE^J J] exp { - \ z* zi/o?} ; (* e C) , (3.3) 

i=i 

the total probability of the 'many' random Gaussian factors immediately gives a severe loss of a total 'relevant 
hitf in phase space with a probability sufficiently close to one. Therefore, it can be advantageous in most cases 
to change to spherical coordinates after a rescaling of the individual random variables Zi —¥ Z{Oi. One thus 
achieves a probability distribution of a radial variable p times a surface area dOaiV,— 1 m the created 2A r 2 -dim. 
Euclidean phase space (with modified scores due to the rescaling with the variances cij) 

N z N z 

II ^T 1 ex P { - \ z * z i/ a i] — ► (* -> z i °i) — > ( II a ^2iV z -i dp p 2 ^" 1 exp{-I p 2 } . (3.4) 
i=l i=l 

If the corresponding score function only depends weakly on the surface term d£l2N z —i, e.g. with lowest order 
spherical harmonics depending only on a 'fevf angular-coordinates, one has obtained an improvement according 
to the principle of 'importance sampling 1 , however even, in the case of a many body physics. 

Further approximation schemes concern the reduction of our chosen, strongly time dependent problem 
of ultra-short laser propagation to the case of a stationary self-energy matrix which is then only subject 

(77 ■ ) (77 ■ ) 

to the difference tj 2 — tj 1 of its two contour time arguments T i2 J2 , (|2,2H2.4b"j) with further combina- 

tion of the corresponding time contour branches r]j 2 , T)j x = db. A Fourier transformation then simplifies 
the whole problem to the case of a self-energy matrix which only entails an individual frequency variable 

, x±) — > S^a^/tiijsi ( w i X2,x\). This simplifies the solution of the saddle point equation 
(I2.14b|) if the physical system is reduced to the case of a stationary time problem. Other approximations and 
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changes of the random variables aE^-j with modification of the probability density ^[aE^cjJ (|3.1cf) . contain- 
ing the determinant DET( [1 — rh] ) + 2 , have similarly to be resolved from the guiding principle of 'importance 
sampling'' in order to acquire a sufficient number of ' relevant hits'' from the random sampling even in the case 
of a many body problem. 
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